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Abstract: In order to compute the entanglement entropy for a given region in a theory with an 
Einstein gravity dual, the Ryu-Takayanagi prescription tells us that we must compute the area of 
an extremal surface anchored to the entangling region. However, if the dual gravity theory receives 
higher-curvature corrections we are compelled to extremize a quantity which is no longer given by 
the area but a higher-derivative functional. Hence, in order to find the extremal surface that yields 
the correct value of the entanglement entropy, we must include an additional boundary condition to 
the problem. We claim that the additional condition can be fixed by demanding that the relationship 
between the bulk depth and the size of the entangling region is the one induced by geodesics, we call 
this the free-kick condition. We implement this prescription in the computation of the entanglement 
entropy of the hairy black hole in new massive gravity, and find a perfect agreement with conformal 
field theory expectations. 
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1 Introduction 

A major force driving much progress in high-energy theoretical physics is the AdS/CFT correspon¬ 
dence or gauge/gravity duality [1]. This duality relates systems that are different in many respects 
by providing an accurate dictionary to relate quantities and phenomena in one system to those in 
the other. More precisely, with the help of the correspondence, we can formulate questions regard¬ 
ing quantum gravity in asymptotically AdS spacetimes as problems in a lower-dimensional gauge 
theory, and vice versa. Moreover, hard problems in one side of the duality are translated into simple 
problems in the other. This feature of the duality allows us to treat analytically questions that were 
hitherto out of the grasp of our tools, e.g. strongly coupled field theories. At the same time, this 
aspect of the correspondence makes it hard to prove its validity in general. In spite of that, the 
AdS/CFT correspondence has been subjected to numerous highly nontrivial tests and it has passed 
them with flying colors. 

Entanglement entropy is a very interesting quantity that can be defined for any quantum system. 
Roughly speaking, the entanglement entropy of a subsystem quantifies the amount of information 
that we would forfeit if we were to lose access to the rest of the system. It must be said, that analytic 
computations of entanglement entropy, in general, can be rather hard. However, there is a great 
body of literature with many interesting results, see, e.g., [2, 3]. In recent years, there has been a 
renewed interest in the study of entanglement entropy. This is partly due to the reformulation of the 
problem, under the light of the AdS/CFT correspondence, by Ryu and Takayanagi (RT) [4]. Their 
proposal has been applied with great success to a wide variety of systems. Moreover, it has been 


- 1 - 




generalized in many different respects; of particular interest for this work has been the extension of 
this prescription to theories whose gravitational dual receives higher curvature contributions [5-9]. 

Three-dimensional gravity has played an important role in the holographic study of entangle¬ 
ment entropy. This is no doubt related to the great deal of analytic understanding that we possess of 
both 3d gravity and 2d conformal field theory (CFT). For example, the entanglement entropy for an 
interval in (l-l-l)-dimensional CFTs has been computed exactly both at zero [10] and at finite tem¬ 
perature [11], these two results have been matched using the RT prescription. Three-dimensional 
higher-derivative gravities have been the subject of ample interest in recent years. In this work we 
focus on one such theory, namely, new massive gravity (NMG) [12]. It is a well-known fact that 
gravitons in three dimensions carry no degrees of freedom; however, if we endow the graviton with 
a mass it will carry two degrees of freedom instead. Thus, we expect a theory such as NMG to have 
attractive new properties to study. Indeed, this theory admits interesting nontrivial solutions and, 
specifically, black holes with an extra gravitational hair parameter [13] . Some properties of these 
black holes have been studied in the context of holographic dualities [14-19]. In the present work, 
we wish to further these studies by investigating in detail the holographic entanglement entropy of 
an asymptotically AdS hairy black hole in NMG. 

Due to the additional hair parameter of these black holes we encounter a puzzle during the 
calculation of holographic entanglement entropy. Indeed, the hair parameter leads to an ambiguity 
in finding the extremal surface used for computing the entanglement entropy holographically, since 
the geodesic is not the entangling curve anymore. The question then is what is the appropriate 
boundary conditions to find the entangling curve? This is the main motivation of this paper. We 
make a proposal for this issue and test it successfully against GFT expectations. We claim that the 
ambiguity can be hxed by demanding that the relationship between the bulk depth and the size of 
the entangling region is the one induced by geodesics, we call this the free-kick condition. 

The plan of this paper is as follows: In Sec. 2 we review massive gravity in (2-|-1)-dimensions 
and then we restrict ourselves to NMG theory. We present the asymptotically AdS hairy black hole 
solution of this theory along its thermodynamical properties. In Sec. 3 we mention the known facts 
about entanglement entropy in (l-l-l)-dimensional GFT and its geometrical interpretation in the 
AdS/GFT context. We then introduce the entropy functional used for computing the entanglement 
entropy in a general higher-derivative gravity theory. We also fix our notations once and for all. In 
Sec. 4 we study the geodesics in an asymptotically AdS hairy black hole background and derive an 
analytic expression for them. In Sec. 5 we use these results and calculate the entanglement entropy 
of the BTZ black hole in NMG. In Sec. 6 we perform an exact analysis of holographic entanglement 
entropy for the asymptotically AdS hairy black hole solution of NMG. We state our proposal to 
impose constraint on the extremal surface used for computing entanglement entropy holographically. 
We then proceed to interpret our results for holographic entanglement entropy from the boundary 
field theory perspective and verify it. We also comment on the closed extremal surfaces in this 
background. In Sec. 7 we summarize our results and discuss the possible future directions. In 
Appendix A, we briefly introduce the special functions that we use in the rest of the paper. 
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2 Review of new massive gravity 


In this work we wish to study the entanglement entropy of a (2+l)-dimensional asymptotically AdS 
hairy black hole. This black hole is a solution of a theory called NMG, which we discuss briefly in 
the present section. Consider the three-dimensional massive gravity described by the action 


S = 


1 


lOvrG 




R-2A + ^CS{T) + ^K 

2/Xcs 


( 2 . 1 ) 


where 


cs(r) = , (2.2) 

K = -Ir\ 

O 

The first term, CS(r), is the gravitational Chern-Simons term of topological massive gravity, while 
K is the higher curvature term introduced in [12] which defines the NMG theory. The field equations 
for (2.1) are 

RfiU 2^9fiu T T 2777,2 = 0 , (2-3) 

where is the Cotton tensor given by 




(2.4) 


and reads 

K^i, = 2\JR^i, - ^V^VuR - hjRgfj,i, + AR^aupR"^^ - - RapR'^'^9^lu + ^R^9iJ.iy, (2-5) 

and fulfills K = 

In the rest of the paper we will consider the parity preserving action, i.e. ^cs co- In the 
following section we will take the value of m such that the NMG admits an asymptotically AdS 
hairy black hole. 


2.1 Asymptotically AdS hairy black hole at L2 = 1/2 

For the special point in the parameter space m?L? = 1/2, L being the AdS radius, the above theory 
(2.1) admits the following black hole solution which obeys weakened AdSa asymptotic boundary 
conditions [13], 

ds^ = —f{r)dt^ + f{r)~^dr^ + r'^dcj)^ , (2.6) 

with 

f{r) = -9 + br+ (2.7) 

Notice that the above black hole solution has vanishing Cotton tensor and hence it satisfies the 
full massive gravity equations of motion (2.3). This black hole is labeled by two parameters; the 
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integration constant n and the gravitational hair parameter h. The mass of the black hole is given 
by 

M = ^. (2.8) 

In the 6 —)• 0 limit the solution is reduced to the BTZ black hole, while in the case of nonvanishing 
h the geometry develops a curvature singularity at r = 0. In the case of 6 > 0 there is a single event 
horizon located at 

r+ = ^ (-bL^ + , (2.9) 

provided the bound fj. > 0 is satisfied. We can associate a Hawking temperature and entropy to it 




( 2 . 10 ) 

( 2 . 11 ) 


Notice that the entropy of the BTZ black hole, i.e., 6 = 0, is twice the one obtained in pure Einstein 
gravity. Henceforth, we refer to (2.6) along with (2.7) as the Oliva, Tempo, Troncoso (OTT) black 
hole. 


3 Holographic entanglement entropy for NMG 

Entanglement is the characteristic trait of quantum mechanics [20] and as such it is important 
to devise concrete methods of quantifying it. Eor example, consider the following simple idea. 
Imagine that the system of interest is in a pure state jT) and suppose that we wish to quantify 
the entanglement between a subsystem A and its complement A. As a first step we could find the 
reduced density matrix obtained by tracing out the degrees of freedom in A 

PA = TTn_^\^){^\. (3.1) 

Once we have constructed this matrix, we notice that if there is any entanglement between the 
degrees of freedom in A and those in A the system appears, to an observer having access only to A, 
to be in a mixed state. If that is the case, then the Von Neumann entropy of pA 

5ee(^) = -TT{pAlogpA ), (3.2) 

is nonvanishing. We refer to this quantity as the entanglement entropy of A. 

In practice it is rather difficult to carry out the procedure outlined above due to technical 
limitations. In spite of this, many general results have been found (e.g., strong subadditivity, area 
law divergence, etc.) and a great deal of progress has been made. Of particular interest to us is the 
fact that in a (l+l)-dimensional GET the entanglement entropy of an interval of length i can be 
computed analytically and it is given by 

S'ee(^) = 2 log , (3.3) 
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Figure 1. The Ryu-Takayanagi prescription; A is the entangling region, and S is the extremal surface in 
the asymptotically AdS background. 


where c is the central charge of the CFT and e an ultraviolet cutoff [10, 11]. Moreover, if we put 
the system at finite temperature and assume periodic boundary conditions, we find 


S'ee(^) 



'A 

vre 


sinh 



(3.4) 


where /3 is the inverse temperature [11]. 

The study of entanglement entropy has been living a period of renaissance during the past few 
years, mainly thanks to the holographic reformulation of the problem due to Ryu and Takayanagi 
[4]. In its original form, the RT proposal claims that for a theory with an Einstein gravity dual the 
computation of the entanglement entropy can be recast as a Plateau problem in an asymptotically 
AdS spacetime. From a practical point of view, in order to compute the entanglement entropy for 
a subsystem A in the boundary theory one needs to extremize the functional 

5ee(^) = ^ (f~^yVh , (3.5) 

where S is a codimension two hypersurface which we demand to be anchored at dA (see Fig. 1), 
and h is the induced metric on S. Hence the entanglement entropy of A is given by 

See{A) = ^ min Area(S). (3.6) 

AG dA=dT. 


Notwithstanding its great efficacy, the prescription (3.6) has two important limitations. First of 
all, it is valid for static spacetimes only, this issue has been addressed by Hubeny, Rangamani and 
Takayanagi by developing a covariant formulation of the RT proposal [21]; in this paper we focus on 
static solutions so we will not discuss the details of this development. On the other hand, experience 
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with black hole entropy [22] teaches us that in the presence of higher curvature corrections, entropy 
computations must be modified; this general wisdom should obviously apply to the RT prescription. 
This question has been discussed recently in [5-9] . In the following we summarize the results in this 
direction that will be relevant for the forthcoming sections. 

Let us consider a general four-derivative gravity action 

5 [R-2A + ciR2 + c2Rf,uR^^ + c^R^upaR^^'^'^] . (3.7) 


For these models, the entanglement entropy is still determined using a codimension two surface 
extending into the bulk as in Fig. 1. However, the shape that S is compelled to take in order to 
yield the correct value of the entanglement entropy is now determined by extremizing the functional 


[5-7] 




1 -|- 2ciR C2 ( i7|| — 2 ^^ ) ^*"3 ( *^11II ~ 


(3.8) 


instead of (3.5). In the following we explain how to obtain the different geometric quantities entering 
into this functional. 

The first step is to hnd a basis for the vector space normal to the surface S. In particular, we 
choose basis vectors and such that 




(3.9) 


where t/q/j is the two-dimensional Minkowski metric. Once we have found these vectors it is possible 
to construct all the ingredients entering Eq. (3.8). The induced metric is simply 


while the relevant contributions from the ambient Riemann curvature read 


(3.10) 


R\\ II = (^(a))'" (^(<5))^^ {^W)Y {n{'l)T ’ (3-11) 

and 

R\\=r^{n^o.)Y{n^p^^R^,. (3.12) 

In turn, the extrinsic curvature is given by 

(^(o))^j, = ) (3.13) 

where V is the covariant derivative wrt The contractions of /C(cj) entering the functional (3.8) 
are 

/C2^r/“/'(/C(„));(/C(^));, (3.14) 

and 

Tr(/C)2 = r?“^(/C(„));(%)),^ (3.15) 


- 6 - 





In the following we shall be concerned with solutions to the NMG theory (2.1). In the notation 
of Eq. (3.7) this theory corresponds to the coefficients 


Cl = 


8m^ ’ 

Therefore the functional (3.8) reduces to 


C2 = 




C3 = 0. 


(3.16) 




r 1 / 

1 . 

3 M 

1 + ^ Rii - -/C" 


\ 

" 2 

4 )\ 


(3.17) 


The computation of entanglement entropy using this functional is the task upon which we embark 
in the following sections. 


4 Geodesics in the OTT geometry 

According to the RT prescription (3.6), the entanglement entropy for an interval A in a (1+1)- 
dimensional GET corresponds to the length of a geodesic extending into the bulk attached to the 
ends of A. In this paper we are interested in computing the entanglement entropy for a system 
dual to the OTT background (2.7); as seen in Sec. 2.1, this background is a solution of a higher 
curvature theory, NMG. Therefore, we do not expect the entanglement entropy and the geodesic 
length to match, instead we expect (3.17) to yield the correct answer. However, the geodesics in 
the OTT background will play an indirect yet important role in the computation of entanglement. 
Therefore, we devote the present section to the study of these curves. 

Geodesics are found by extremizing the functional 

I{A)= [ d^-^yVh, (4.1) 

Jt.a 

where A is a spacelike interval at asymptotic infinity and has its endpoints fixed at dA. We 
could either parametrize these curves by the angular or by the radial AdS coordinates, hereafter we 
refer to the former as the boundary parametrization and to the latter as the bulk parametrization. 
Ghoosing the bulk parametrization, we find that the induced metric is given by 

h = - /(r)"^, (4.2) 


where /(r) is given by Eq. (2.7). Notice that the functional (4.1) does not depend explicitly on (f> 
but only on its derivatives. Setting L = 1, we find 


0'(r) 


1 

[r — r+)(r + r_|_ + b){ar‘^ — 1) ’ 


(4.3) 


where a is an integration constant and r_|_ is the horizon of the black hole (2.9). Therefore, the 
profile of a geodesic reaching down to a radius r* > r+ is given by 


(/>(r) 



df . 


(4.4) 
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The integration constant a can be determined by imposing the boundary condition 


dr 

d(j) 


= 0 , 


(4.5) 


which implies that a = r^‘‘. Hence, the bulk-parametrized geodesics read 


4>{r) = j 
J r-. 


= I dr 

r* r^y{r - r+){f + r++b){r‘^ - r'^) 


(4.6) 


Performing this integration explicitly, we find 

2 r^F{z\r]) + (r+ — r*) n(n; z\r]) 


(j){r) = 


v^(r+ + r*)(6 + r+ -b r*) 


(4.7) 


where F{z\ri) and n(n; 2 ;|r^) are incomplete elliptic integrals of the first and the third kind, respec¬ 
tively and 


2r4 


n = 


r+ -I- r* 


z = arcsm 


(r - r*)(r+ -h r*) 


V = 


2r*(r — r+) 

2r*(6 -I- 2r+) 

(r+ r*)(6 r+ -b r*) ' 


(4.8a) 

(4.8b) 

(4.8c) 


More details about these integrals can be found in Appendix A. 

Notice that in the bulk parametrization the depth r* reached by the geodesic is the tunable 
parameter (see Fig. 3), from the boundary point of view; however, one would wish to specify the 
size of the entangling region (p instead. It is clear that these two quantities are related by 


'^geodesic(^*) 


(4.9) 


To simplify notation, we drop the subscript “geodesic” when no confusion arises, and unless oth¬ 
erwise stated. Regretfully, it is quite complicated to solve this equation in general; however, it is 
possible to do so by setting 6 = 0, in this case the OTT solution becomes the BTZ black hole and 
(4.7) reduces to (see Fig. 2 to compare the OTT and BTZ geodesics) 


4>{r) = — arccosh 

r+ 


9 9 

rt — r± 


(4.10) 


Here, it is straightforward to go from one parametrization to the other and the boundary parametriza¬ 
tion is simply 

cosh^ {r^4>) \ ^ 


1 - 


cosh^ (r+0o)/ 


r((/>) = r+ 


(4.11) 

















Figure 2. The comparison between geodesics in the BTZ (dashed) and OTT (solid purple) geometries; the 
radial coordinate is compactified using the map r —> arctan(r). The outer (solid black) circle represents 
asymptotic infinity, and the inner (dashed) circle the horizon radius. 


where 



(4.12) 


Using Eq. (4.9) we observe that 4> = (po- If 6 yf: 0 this relationship is more complicated but the 
relevant information can be extracted systematically, we will come back to this issue in Sec. 6. 


5 Entanglement entropy of the BTZ black hole in NMG 


Before attempting to compute the entanglement entropy of the OTT black hole, we set 6 = 0 and 
perform this computation for the BTZ black hole. The first step is to find the geometric quantities 
entering Eq. (3.17), the induced metric h and i2|| are given by 


h = 


1 




Rii = — 


while the contraction (3.14) of the extrinsic curvature reads 


L2 ’ 


(5.1) 


K^ir) = 


2L?‘r + 2r‘^ + [6L^r‘^ — 4^L^) 


4L2 [r2(r2 - ^7.2) </.'(r)2 + L2]' 


(5.2) 


Plugging these into (3.17) and carrying out the variation we find a complicated equation of motion 
for 4‘{r). Interestingly enough, the solutions of the geodesic equation (4.10) solve also these more 
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Figure 3. Conformal diagram of entangling curves for different depths in the BTZ geometry. 


convoluted equations of motion as pointed out in [23] ^ . 

Notice that although the relevant extremal curves are the same as in Einstein gravity it is not 
their length that we must consider, instead we ought to insert (4.10) into (3.17) and compute the 
integral. By doing so, we find that the entanglement entropy of the BTZ black hole in NMG is 
given by, setting T = 1, 


See = 


1 


dr 


^ (r2 - rl) (r2 - rl) 

llog(^r2-q + \A^) 


(5.3) 


where we introduced an ultraviolet cutoff 1. Finally, using (4.12), we can rewrite this expression 
in the boundary parametrization. We replace r* using 


= 

I ^ 


1 — cosh ^ ( r_|_ 


and obtain 


5'ee = ^ log 


—^ sinh ( r_(_ 
r+ 


Finally, introducing the quantities 


f3 = — 
r+ 


r. 


£ = 24,, 


(5.4) 


(5.5) 


(5.6) 


^For the BTZ geometry, the higher-derivative terms contribution is topological and can be written as a total 
derivative. Therefore, the holographic entanglement entropy is still determined by an extremal length curve in the 
bulk. 
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we find 


5'ee = ^ log 


/3 


sinh 


7re 


(?) 


(5.7) 


If we compare this result with (3.4) we find that the expressions match if 


3 


(5.8) 


The reader might find this puzzling since this c seems to be twice the value of the central charge 
computed by Brown and Henneaux [24]. However, this contradiction is only apparent since the 
central charge for NMG is given by [16, 17] 


c = 


3 




(5.9) 


which at the special point = 1/2 is precisely (5.8). Given that one can always find a BTZ black 
hole solution in NMG regardless of the point in parameter space, the computation above can be 
replicated for arbitrary m; doing so, we find 


See 


1 



'A 

vre 


sinh 



(5.10) 


as expected. It is pleasant to see how (3.17) captures this shift of the central charge in a nontrivial 
manner. 


6 Entanglement entropy of the OTT black hole 


The time for computing the entanglement entropy of the OTT black hole has arrived. This com¬ 
putation will exhibit all the intricacies of calculating entanglement entropy in a higher curvature 
theory, yet it is simple enough so that analytic expressions can be found. The entropy corresponding 
to closed entangling curves was computed in [23]. Below, we focus on the case where these curves 
are anchored at infinity. We shall see that to deal with this case it is necessary to introduce a 
new prescription to determine the boundary conditions. Once again, we kick off by computing the 
geometric quantities entering (3.17). The induced metric h and i?|| are given by 


h = 


+ rV(r)^ 


^ (^+ A) “ ^ 


i?ll = — 


bL^ 


2r [r^(/>'(r)^ (^L^r —-|-r^) + L^] r ’ 


( 6 . 1 ) 


while the extrinsic curvature contraction (3.14) reads 
}C^(r) = ^ 


4L^ [r20'(r)^ {hL'^r — + r^) + ' 

„ 2^2 


2L^r(/)"(r) {pL^r — 


-|- 2r^(/'(r)^ {pL^r — iiL? -[- r^) + p'{r) {pbL'^r — + 6L^r^) 


( 6 . 2 ) 
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For the present case, the Euler-Lagrange equation we need to solve is given by 


(f f 5C 


d f 5C 




5C 


5(f) 


(6.3) 


where C follows from (3.17), (6.1) and (6.2). Moreover, since this Lagrangian does not depend 
explicitly on (f>{r), we have 

d f 5C\ 5C ^ 

V^y ■ 

where k is an integration constant. Putting everything together, the equation of motion for (j){r) is 
given by the following complicated expression, 

1 




X I - L'^r‘^(t)'{rf (dl^^L"^ + 805LV + - 4/xLV [TihL^ + 20r) + 56/r^L^ 

+ 4:r^(f>'{ry (/iE^ + r^) {hL^r — fiL^ + + Ar^(j)'{r)^ (^—bL^r + QfiL^ + 2r^) 


Lr {hL?' + r) — 


1 2 


+ 4L®r 


2r(f)^^\r) (bL^r — + r^) 


+ (f)"{r) {flbL^r — AgiL^ + +AL^(j)'{r) — f)"{r)"^ {pL^r — gL? + 

+ 2bL^r + AgL'^ + + SL^r^cf'{r)‘^ {bL^r — fiL^ + r^) r(f>^^\r) {bL^r — + r^) 

+ (j)”{r) iy—AblF'r + 3/iL^ — 5r^) > = /c. 


(6.5) 


In the following, we proceed to solve Eq. (6.5) to linear order in b. Therefore, we write 


0(r) = + b(f)^^\r) + O (6^) , 

k = +bk^^^ + 0{b‘^) , 


( 6 . 6 ) 

(6.7) 


and use (2.9) to express g, in terms of the horizon of OTT black hole and b. Then we expand (6.5) 
around 6 = 0. Obviously, the zeroth order contribution reduces to the NMG-BTZ black hole, hence, 
<^( 0 )(r) is given by the geodesic prohle (4.10) and k^^^ = 2r*. Nontrivial deviations from the geodesic 
path appear already to linear order in 6 , setting L = 1 , we find 


1 ^( 1 ) (r) = 



r — r+ (r + r+) — r* (rj + rr+) 


r + r+ 


2rl {r\ - rf) 


+ Cl 




(/>(°)(r) + C 3 . 


( 6 . 8 ) 


Observe that ( 6 . 8 ) has four unknown constants ci, C 2 , C 3 , and k^^'^ that we ought to fix. Clearly, 
this is a consequence of the fact that the equations of motion (6.3) is a fourth order differential 
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equation, but the pressing question is: How are we going to fix these extra constants? 

For a moment let us go back to the bulk parametrization. Obviously we are also expected to 
fix four independent constants there. Three of them are rather easy to fix by setting 




= C* , 


. =0’ 

r'"(0) 

</)=0 


<t>=0 



= 0 , 


(6.9) 


we are already familiar with the first two conditions while the third one can be motivated by 
symmetry. The missing parameter to fix is the initial acceleration or the concavity of the entangling 
curve. In Fig. 4 the different solid curves represent solutions complying with (6.9) for different 
choices of r"(0). Hence, there seems to be an unfixed parameter in the game, so the natural 
question now is whether there is any preferred initial acceleration. 

It is at this point that the discussion in Sec. 4 becomes relevant. We propose that this freedom 
can be fixed by demanding that the curve used for computing the entanglement entropy holographi¬ 
cally must not intersect the shadow enclosed by the geodesic (4.7) (this is the gray region in Fig. 4). 
More explicitly, we claim that the acceleration at (/> = 0 has to be chosen such that 


— <^geodesic(^*) • 


( 6 . 10 ) 


This corresponds to the dashed line in Fig. 4. Hereafter, we refer to (6.10) as the free-kick condition 
since it is analogous to the problem that a football player has to solve when executing a free-kick. 
Namely, given an object at rest at a fixed position, what is the acceleration necessary to hit a 
predetermined target. In this analogy, the player is placed at the tip of the entangling curve, while 
the target is given by the asymptotic value of the geodesic at spatial inhnity. Back in the boundary 
parametrization, the boundary conditions (6.9) read 


1 


= 0 , 


= 0 , 


0 '(r) 

‘i(j)"{r)'^ — (f)'" {r) 

(/)'(r)^ 


r=r* 


= 0 . 


From (6.8) we see that the conditions (6.10) and (6.11) are satisfied by setting 

.2 _ ^( 1 ) 2 


Cl = 


2rl^rl-rl 

C2 = C3 = 0 , 

^( 1 ) _ (c-H - r^Y 


r+ 3r+(r++r*)3/2 


2 r* 


+ 


r. , 3 1 3 5 
' 2 ’ 2 ’ 2 ’ 

2 


(r+ — rY) r^F{w\p) + (r+ -|- r*) E{w\p) 


ri (r+ -I- r*) 


(6.11a) 

(6.11b) 

(6.11c) 


( 6 . 12 ) 
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Figure 4. Conformal diagram of extremal curves in the OTT geometry for different concavities (solid purple), 
geodesic curve shadow (gray), and the free-kick extremal curve (dashed). 


where E is the incomplete elliptic integral of the second kind (A.lb), Fi is the Appell hypergeometric 
function (A.2), and 


w = arcsin 


/ /r+ -h r* 

VV 2r* 



P = 


4r+r* 
(r+ -h r*)^ 


Additionally, we find that to linear order in b the acceleration is given by 


r"i<P) 



(j)"{r) 


= [rl - r+) -h ^ - 2r+ -h 6. 

r=r* 


(6.13) 


(6.14) 


Finally, we insert (4.10), (6.8), and (6.12) into (3.17) and find the entanglement entropy of the 
OTT black hole to linear order in b 


SEEir+,r^) = 


2Grl 


b\/r'^ ~ r+ (^- 1 - ~ + r (rf — -|- r+rf — 


r (r + r+) (r^ - 4) 


s/r'‘-r‘i^rl-rl V / ^ 


+ 0(6^). (6.15) 


As can be seen from the above, the holographic entanglement entropy depends on the constant 
at first order in b, which cannot be fixed unless extra boundary conditions are provided. Notice 
that, for the BTZ geometry, 6 = 0, the holographic entanglement entropy does not depend on k^^\ 
and therefore we do not need extra boundary conditions; this is consistent with the discussion in 
Sec. 5. This expression seems to have a quite complicated dependence on the bulk parameters r+ 
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and r*; we would like to know how (6.15) reads in terms of the boundary parameters fd and (p. We 
test whether this expression matches the CFT expectation 


I log 


^ 

— smh — 
7re V 


(6.16) 


with /3 the inverse temperature of (2.10) and c the central charge (5.8) corresponds to the sought 
after S'ee(/35 ^)- We would like to emphasize that it happens since we have fixed the constant 
using our proposed free-kick condition. It is clear that other curves do not fit the CFT expression 
for the entanglement entropy of an interval at finite temperature, since the ordinary differential 
equation (ODE) is fourth order and we used four boundary conditions to fix these freedoms for 
finding the entangling curve. We proceed to verify this claim. First, we set (. = 2(j) m. (6.16) and 
use (6.10), i.e., 

2 r^F{z\ri) + (r+ — r*) n(n; z\ri) 


(pin) = 


(6.17) 


'C+ ■\/(r+ -|-r*)(6-Fr+ + r*) 

where 5 = limj._>.oo -Z; and z is defined in (4.8b). Then we express /? in terms of the horizon r_|_, using 
(2.9) and (2.10) we find 


/3 = 


dvr 


(6.18) 


6^ + 2r_|_ 

Inserting the above expressions into Eq. (6.16) and expanding around 6 = 0 we find that it matches 
(6.15) exactly at first order in b, this gives strong support to our claim. 

Einally, we verify our claim numerically by adapting the algorithm discussed in [25] to the 
present case. In order to get useful numerical results, we must compute the renormalized version of 
entanglement entropy [26, 27], which for a (l+l)-dimensional CFT is given by 


^EE = f’v^'S'EE • 

Ol 

The renormalized entanglement entropy corresponding to (6.16) reads 


c£ ( 


(6.19) 


( 6 . 20 ) 


Running a code that implements the free-kick condition alongside with the algorithm introduced in 
[25] we find that the 5 ee resulting from solving Eq. (6.5) matches (6.20) within numerical error. 

To finish this section, we briefly comment on the closed extremal surfaces in the OTT back¬ 
ground. The black hole bifurcation surface r_|_ is a closed extremal surface of the entanglement 
entropy functional (3.17), and the corresponding entanglement entropy is the Wald entropy of the 
black hole 




EE 




( 6 . 21 ) 


If we set cj)' = (j)" = (f)'" = (()'"' = 0, we find another closed extremal surface given by 


ra = , 


( 6 . 22 ) 
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whose corresponding EE is simply 

• (6.23) 

Notice that r+ > Tq. In the absence of the gravitational hair parameter b, i.e., the BTZ black hole, 
the bound is saturated and is exactly the Wald entropy of the BTZ black hole. 

7 Conclusions and outlook 

Since pure Einstein gravity in (2+1) dimensions admits few solutions, one is compelled to consider 
more exotic actions like NMG. The NMG theory admits a rich and interesting catalog of solu¬ 
tions such as black holes, wormholes, solitons, and kinks. Hairy black holes are a characteristic 
of NMG and one would like to have a holographic interpretation for them. The OTT black hole 
that we considered in this paper has the intriguing feature that geodesic curves do not correspond 
to the extremal surface used for computing the entanglement entropy holographically. So, this 
work has established that if one would like to calculate the holographic entanglement entropy in 
higher-derivative gravity theories one should first impose additional boundary conditions to hnd the 
extremal surface. We demonstrated that this can be done by demanding that the entangling curve 
must not intersect the region enclosed by the geodesic. 

We then computed the holographic entanglement entropy of the OTT black hole. We conhrmed 
that the extremal surface which is hxed by proposing our boundary conditions yields the known 
result for entanglement entropy in the (l+l)-dimensional GET at finite temperature. 

In the recent paper [28], the authors considered AdS wave solutions of NMG which admit 
logarithmic modes in their Eefferman-Graham expansion. These solutions are believed to be dual 
to logarithmic CETs. They calculated the holographic entanglement entropy for these backgrounds 
and showed that the entanglement entropy has a new divergent term. Beside the aforementioned, the 
NMG theory at the chiral point, where central charges vanish, admits a logarithmic deformation of 
extremal BTZ as an exact solution. The entanglement entropy for a GET dual to the extremal BTZ 
black hole was computed recently [29] . It is an interesting future problem to study the implications 
of the free-kick boundary conditions to address the logarithmic deformation of extremal BTZ in 
boundary conformal field theory. 

In this paper, to avoid even more cumbersome computations we focused on the static case; 
however, it would be plausible to extend our current study to the rotating hairy black holes. Due 
to the gravitational hair parameter the asymptotically AdS rotating hairy black holes have two 
different kind of extremal limits. In one of them, both the temperature and the entropy vanish 
and, therefore, it defines a Nernst solution. It will be interesting to investigate the behavior of 
holographic entanglement entropy for these geometries. 

Lastly, we would like to comment on some recent work which showed that one can put constraint 
on the extremal surfaces by imposing causality [23]. We proposed new boundary conditions to 
completely fix the extremal surface used for computing holographic entanglement entropy in higher- 
derivative gravity theories. One of the intriguing future directions is to try to understand the relation 
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between the argument based on CFT causality [30] and our work. These are the avenues that we 
will try to explore in the near future. 


A Elliptic integrals 


In this section we will shortly review the definitions of the special functions which we used in this 
paper. The incomplete elliptic integrals can be written as 

^ r dt 
F{z\v)= /■ . 

JO yl —r/sm (t) 

E{z\r]) = J — r] sin^ {t) dt , 

TT/ , ^ r dt 

{n,z\r])-J^ -^=====j——^^ 

The parameter r] is called the modulus of the elliptic integral and z is the amplitude angle. They 
range from 0 < z < ^ and 0 < r] < 1. 

The Appel hypergeometric function Fi is defined by the double series 


(A.la) 
(A.lb) 
(A.lc) 


Fi (a; bi,b 2 ; c; zi, ^ 2 ) = 


{bi)k {b2)i ia)k+i 

k\l\{c)k+i 


\zi\ < 1 A IZ 2 1 < 1 , 


(A.2) 


k=0 1=0 

where the Pochhammer symbol {q)n denotes the rising factorial, 


{q)n = q{q+l)...{q + n-l). 


(A.3) 


Acknowledgments 

We are grateful to Shajid Haque for commenting on a draft of this paper. We would like to thank 
Mohsen Alishahiha, Arpan Bhattacharyya, Pawel Caputa, Amin Faraji Astaneh, Mario Flory, Kevin 
Goldstein, Charlotte Sleight, and Amir Hadi Ziaie for illuminating conversations. The research of 
A.V.O. is supported by the University Research Council of the University of the Witwatersrand. 
The work of S.M.H. is supported in part by INFN. S.M.H. and A.V.O. wish to thank the hospitality 
of the ICTP in Trieste where part of this research was conducted. A.V.O. is grateful to the theory 
group at Universita di Milano-Bicocca for kindly hosting him during the concluding stages of this 
work. 

References 

[1] J. M. Maldacena, “The large N limit of superconformal field theories and supergravity,” Int. J. Theor. 
Phys. 38 (1999) 1113 [Adv. Theor. Math. Phys. 2 (1998) 231] [hep-th/9711200]. 

[2] H. Casini and M. Huerta, “Entanglement entropy in free quantum field theory,” J. Phys. A 42 (2009) 
504007 [arXiv:0905.2562 [hep-th]]. 


-- 17 - 








[3] P. Calabrese and J. Cardy, “Entanglement entropy and conformal field theory,” J. Phys. A 42 (2009) 
504005 [arXiv:0905.4013 [cond-mat.stat-mech]]. 

[4] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from AdS/CFT,” 

Phys.Rev.Lett. 96 (2006) 181602 [hep-th/0603001]. 

[5] J. Camps, “Generalized entropy and higher derivative Gravity,” JHEP 1403, 070 (2014) 

[arXiv: 1310.6659 [hep-th]]. 

[6] D. V. Fursaev, A. Patrushev and S. N. Solodukhin, “Distributional Geometry of Squashed Gones,” 
Phys. Rev. D 88 (2013) 4, 044054 [arXiv:1306.4000 [hep-th]]. 

[7] X. Dong, “Holographic Entanglement Entropy for General Higher Derivative Gravity,” JHEP 1401, 
044 (2014) [arXiv:1310.5713 [hep-th], arXiv:1310.5713]. 

[8] A. Bhattacharyya, M. Sharma and A. Sinha, “On generalized gravitational entropy, squashed cones 
and holography,” JHEP 1401, 021 (2014) [arXiv:1308.5748 [hep-th]]. 

[9] S. N. Solodukhin, “Entanglement entropy, conformal invariance and extrinsic geometry,” Phys. Lett. B 
665 (2008) 305 [arXiv:0802.3117 [hep-th]]. 

[10] G. Holzhey, F. Larsen and F. Wilczek, “Geometric and renormalized entropy in conformal field 
theory,” Nucl. Phys. B 424 (1994) 443 [hep-th/9403108]. 

[11] P. Calabrese and J. L. Cardy, “Entanglement entropy and quantum field theory,” J.Stat.Mech. 0406 
(2004) P06002 [hep-th/0405152]. 

[12] E. A. Bergshoeff, O. Hohm and P. K. Townsend, “Massive Cravity in Three Dimensions,” Phys. Rev. 
Lett. 102, 201301 (2009) [arXiv:0901.1766 [hep-th]]. 

[13] J. Oliva, D. Tempo and R. Troncoso, “Three-dimensional black holes, gravitational solitons, kinks and 
wormholes for BHT massive gravity,” JHEP 0907, 011 (2009) [arXiv:0905.1545 [hep-th]]. 

[14] R. Fareghbal and S. M. Hosseini, “Holography of 3D Asymptotically Flat Black Holes,” Phys. Rev. D 
91 (2015) 8, 084025 [arXiv:1412.2569 [hep-th]]. 

[15] S. de Buyl, S. Detournay, G. Giribet and G. S. Ng, “Baby de Sitter black holes and dS 3 /CFT 2 ,” 

JHEP 1402, 020 (2014) [arXiv:1308.5569 [hep-th]]. 

[16] G. Giribet, J. Oliva, D. Tempo and R. Troncoso, “Microscopic entropy of the three-dimensional 
rotating black hole of BHT massive gravity,” Phys. Rev. D 80, 124046 (2009) [arXiv:0909.2564 
[hep-th]]. 

[17] E. A. Bergshoeff, O. Hohm and P. K. Townsend, “More on Massive 3D Gravity,” Phys. Rev. D 79 
(2009) 124042 [arXiv:0905.1259 [hep-th]]. 

[18] G. Giribet and M. Leston, “Boundary stress tensor and counterterms for weakened AdSa asymptotic 
in New Massive Gravity,” JHEP 1009, 070 (2010) [arXiv:1006.3349 [hep-th]]. 

[19] Y. Kwon, S. Nam, J. D. Park and S. H. Yi, “Holographic Renormalization and Stress Tensors in New 
Massive Gravity,” JHEP 1111, 029 (2011) [arXiv: 1106.4609 [hep-th]]. 

[20] E. Schrodinger, E, “Discussion of Probability Relations Between Separated Systems,” Proceedings of 
the Gambridge Philosophical Society, 31: 555563; 32 (1936): 446451 

[21] V. E. Hubeny, M. Rangamani and T. Takayanagi, “A Govariant holographic entanglement entropy 
proposal,” JHEP 0707 (2007) 062 [arXiv:0705.0016 [hep-th]]. 


- 18 - 



[22] R. M. Wald, “Black hole entropy is the Noether charge,” Phys. Rev. D 48, 3427 (1993) 
[gr-qc/9307038]. 

[23] J. Erdmenger, M. Flory and C. Sleight, “Conditions on holographic entangling surfaces in higher 
curvature gravity,” JHEP 1406, 104 (2014) [arXiv: 1401.5075 [hep-th]]. 

[24] J. D. Brown and M. Henneaux, “Central Charges in the Canonical Realization of Asymptotic 
Symmetries: An Example from Three-Dimensional Gravity,” Commun. Math. Phys. 104 (1986) 207. 

[25] A. Bhattacharyya, S. Shajidul Haque and A. Veliz-Osorio, “Renormalized Entanglement Entropy for 
BPS Black Branes,” Phys. Rev. D 91, no. 4, 045026 (2015) [arXiv:1412.2568 [hep-th]]. 

[26] H. Casini and M. Huerta, “A Finite entanglement entropy and the c-theorem,” Phys. Lett. B 600 
(2004) 142 [hep-th/0405111]. 

[27] H. Liu and M. Mezei, “Probing renormalization group flows using entanglement entropy,” JHEP 1401 
(2014) 098 [arXiv: 1309.6935 [hep-th], arXiv:1309.6935]. 

[28] M. Alishahiha, A. F. Astaneh and M. R. M. Mozaffar, “Entanglement Entropy for Logarithmic 
Conformal Field Theory,” Phys. Rev. D 89, no. 6, 065023 (2014) [arXiv:1310.4294 [hep-th]]. 

[29] P. Caputa, V. Jejjala and H. Soltanpanahi, “Entanglement entropy of extremal BTZ black holes,” 
Phys. Rev. D 89, no. 4, 046006 (2014) [arXiv:1309.7852 [hep-th]]. 

[30] M. Headrick, V. E. Hubeny, A. Lawrence and M. Rangamani, “Causality & holographic entanglement 
entropy,” JHEP 1412, 162 (2014) [arXiv:1408.6300 [hep-th]]. 


19 - 



